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ABSTRACT 

Aims. We study the two-dimensional distribution and kinematics of the stellar and gaseous components in the centre of the prototype 
Wolf-Rayet blue compact dwarf galaxy. He 2-10. The aim is to compare the kinematics of gas and stars in order to determine whether 
they are consistent with one another, or if stars and gas can be decoupled due to gravitational perturbations and feedback from star 
formation. 

Methods. We have used the integral field unit ARGUS, part of FLAMES on the European Southern Observatory's Very Large 
Telescope, to target the Can /1/18498, 8542, 8662 A triplet in the central 300 x 480 parsecs of He 2-10. The selected wavelength 
regime includes several prominent spectral features, including the Paschen series and the [S iii] emission-line, which we have used to 
derive the kinematics of the ionised interstellar medium. 

Results. We find no systematic trend in the velocities of the stars over the observed field of view and conclude that the stellar 
kinematics is governed by random motions. This is in contrast to the motions the ionised interstellar medium, where we find spatial 
velocity variations up to 60 km s~^ . Our gas velocity field is consistent with previous studies of both the molecular gas and the 
feedback-driven outflow in He 2-10. We interpret the kinematic decoupling between the stars and the gas as He 2-10 being in the 
process of transformation to a dwarf elliptical galaxy. 
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. 1. Introduction 

Blue compact galaxies (BCGs) are important laboratories for 
I testing existing theories for dwarf galaxy formation and evolu- 
tion. This is because many of them form stars at rates that are 
unsustainable over long periods and because of their low inter- 
] stellar medium (ISM) metallicities. The large amounts of young 
stars in BCGs give rise to a blue continuum and strong emission- 
lines, dominating the optical light. Their high star formation 
rates, high gas content, chaotic morphology and low metallicities 
make them o ften being r egarded as local counterparts of young 
galaxies (e.g. jKunth & Qs tlin 2000). However, detection of faint 
old stellar populations suggests that these syste ms are intrinsi- 
cally old, even in the most youthful cases (Qstl in & Mouhcin^ 
^2005). The puzzle of what triggers starbursts in BCGs remains 
unsolved, although merger events are likely to play an important 
role. 

Previously, E abry-Perot interferometric studies of the Ha 
line (lOstlin et al . 1999, 2001) have revealed disturbed kinemat- 
ics, indicating dwarf galaxy mergers as a likely explanation for 
the formation of BCGs. Despite local turbulence, the Ha kine- 
matics nevertheless suggest overall gravitational support. Due to 
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the ubiquitous presence of shocks, supernova winds and heating 
by massive stars, it is unclear if the ISM follows the motions of 
the stars, which are bound to follow the gravitational potential. 
Both mass estimates from gas kinematics and evolutionary sce- 
narios based on the assumption that BCGs are self-gravitating 
may therefore be uncertain. This is of importance not only for the 
local universe, but also at higher redshift, where strong emission- 
lines are often the only usable kinematical tracers. Thus, mass 
estimates of high-redshift star-forming galaxies may change sig- 
nificantly if the stars and gas display diflferent kinematics. 

Stellar motions are usually studied using optical absorp- 
tion lines. However, in BCGs, where metallicities are low and 
star formation is strong, the absorption lines become easily di- 
luted by the strong blue continuum of hot stars and the nebular 
emission. In the near-infrared, the Can AAM9^, 8542, 8662 A 
triplet (hereafter Ca triplet) can be used t o derive reliable stel- 
lar kinematic information ( Dresslej|l984^ . In a previous study 
using Very Large Telescope (VLT) long-slit spectroscopy, we 
confirmed on ESO 400-G43 that the Ca triplet can be u sed to 
probe the stellar kinematics in BCGs ( Qstlin et al.ll2QQ4l) . Qur 
study showed that the stellar component did not show the super- 
keplerian decline observed for the ionised gas. Thus, stars can 
be kinematically decoupled from the gas. With the help of inte- 
gral field units (lEUs), we can compare the gaseous and stellar 
motions by spectroscopy in two spatial dimensions. 
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30 arcsec =1.3 kpc 



Fig.l. HST/WFPC2 F814W archival image of He 2-10 with 
an ov erlaid CO (thick white) contour from (Kobulnick v et alJ 
The thin black lines outline the positions of bubbles ob- 
served in Hot emission (iMendez et al .11 19991) . The inset shows 
the central region at larger scale together with the footprint of 
our ARGUS observations and regions A and B, as discussed in 
the text. 



Here, we present results from IFU observations of the central 
300x480 parsecs in He 2-10, targeting the Ca triplet. He 2-10 is 
a well-known BCG with an absolute 5-band magnitude of -19 
and two prominent star-forming regions (see Fig[T]). The super 
star clusters (SSCs) in the cent ral region (region A) have an age 
of 4-5 Myr ( Chandar et al.ll200 3) while those in region B, lying 
east of the centre, are older. It has been argued by Mendez gt^ 
(Il999h that strong feedback from star formation drives a bipolar 
outflow in He 2-10. The systemic velocity of 870 kms"^ is too 
low to serve as a good distance measure because it cannot be 
assumed to belong to the Hubble flow, due to the influence of 
the Virgo cluster Here, w e adopt the distance of 9 Mpc used by 
iKobulnicky et al.l (Il995|) . corresponding to a physical scale of 
46.3 pc arcsec"^. 



2. Observations and Data Reduction 

The observations were carried out under good weather con- 
ditions during two nights in November 2004 with FLAMES 
at European Southern Observatory's VLT (VLT/UT2) and 
its IFU, ARGUS, c onnected to the spectrograph GIRAFFE 
(iPasquini et al.ll2002h . The coarse (1:1) scale of ARGUS gives 
a spatial resolution of 0^.^52 pix 30 pc/pix and the 22 x 14 
pixels result in a field-of-view of 1 V!4 x 7! 3 placed with 40° po- 
sition angle on region A in He 2-10 (see Fig[T]). The light from 
each resolution element is fibre-fed to the spectrograph where 
spectra with R ^ 10400 ranging from 8200 A to 9380 A were 
measured. The spectra therefore cover the Ca triplet, O i /18446, 




8600 



8700 



Fig. 2. Example of two spectra from single fibres with diff'erent 
S/N. Dotted line: original spectrum. Solid line: after subtraction 
of the Paschen lines. 



[Chi] ^8579, [Feiii] ^8617, [Siii] ^906^1 and the Paschen 9 to 
Paschen 19 emission-lines of H i. 

To save overhead time, we used fast target acquisition, 
thereby trusting the pointing accuracy of the telescope. Four 
frames were obtained with a total exposure time of 3.2 hours. 
The data were reduced u sing the software and methods by 
iPiskunov & Valentil (l2002h . with reduction steps including stan- 
dard techniques for bias subtraction, flat-fielding, wavelength 
calibration and optimised extraction of each spectrum. The back- 
ground emission from the sky was then subtracted using the si- 
multaneously observed spectra from the 15 sky fibres lo cated at 
a large distance around the target (see lKaufer et al.ll2003[ for de- 
tails). 

Since all three lines of the Ca triplet have Paschen emis- 
sion lines superimposed on their red wings, it is crucial to accu- 
rately subtract the emission features before deriving the stellar 
kinematics. By fitting several isolated Paschen lines simultane- 
ously and using the theoretical relative line- strengths, we sub- 
tracted a model spectrum for Paschen emission for each pixel. 
The model spectrum was constructed by fitting the uncontami- 
nated lines of Paschen 9, 10, 14, 17, 18, and 19. Manual adjust- 
ments were made in cases where the signal was too low to get a 
reliable model for the Paschen lines. Fig. [2] shows two examples 
of subtracted and unsubtracted spectra in the region around the 
Ca triplet. Both emission and absorption lines are indicated. 

We parametrise the stellar line-of- sight velocity distribution 
by a Gaussia n profile using the penali zed pixel fitting algorithm 
developed by' Cappellari & EmsellemI (12004) toget her with a se- 
lection of template stars from the stellar library by Cenarro et al.l 
(2001). Using this technique, we derived the kinematics using 
the best linear combination of 1 8 carefully chosen template stars 
covering M6 V with Teff =3721 to 09 V with Teff =36300 K, and 
metallicities ranging between -2.25 and 0.13 solar. The stellar 
types used were KO V, A8 Vn, KO III, K3 lib, K4 lab, M3 II-III, 
K3 lab, 09 V, B5, K3 III, A4, KO V, M6 V, M2 V, G5 IIIw, 
K7 V, FO and M4 V. Thus our derived kinematics is not af- 
fected by template mismatching eff'ects, since a diff'erent com- 
bination of stellar template spectra is compiled for each individ- 
ual spectrum. To match the template stellar library, we degraded 
the spectral resolution of our data from 0.2 A/pix to 0.85 A/pix. 
This simultaneously increased the signal in the fainter pixels, 
while still resolving the Ca triplet well enough to derive the 
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dRalchenko et al.|[2007h . corresponding to 16.5 km s 
paper, we adopt 9069.0 A. 



. Throughout this 
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kinematics accurately. However, for fitting the emission lines, 
we used the full-resolution spectra. 

Although the pixel fitting routine delivers formal errors, 
these could be underestimated since this formalism suppresses 
the noise in the observations. We used Monte Carlo simulations 
to calculate more realistic error estimates. For each pixel, we 
derived the kinematics for 250 realisations, i.e., Gaussian ran- 
domised error- spectra added to the galaxy spectra. We found that 
the average error over all pixels is 6 km s"^ for the velocities, and 
9 km s"^ for the stellar velocity dispersions. We estimate that the 
emission line kinematics can be derived within the same level of 
uncertainty or better. 

In order to test the robustness of the derived stellar kinemat- 
ics, we were able to reproduce the stellar velocity field inde- 
pendently by instead cross-correlating with three template stars 
of diff'erent type (KOII, KOIII, G8Iab) that were observed in the 
same nights with the same instrumental setup. To quantify the 
systematic error from the subtraction of the Paschen emission, 
we varied the strength of the subtraction from twice the adopted 
amount to no subtraction at all. The goodness of the subsequent 
fitting confirmed that the subtraction is optimal within 20%. This 
corresponds to a possible systematic error of 10 and 15 kms~^ 
for the stellar velocities and line widths, respectively. However, 
this is probably quite conservative since there is no indication 
that this error acts always towards under- or oversubtraction and 
the spread in stellar velocity values over the whole field is only 
8kms"^ 

3. Results and Discussion 

3.1. Stars 

The Ca triplet is one of the most prominent near-infrared fea- 
tures of cool stars. Although it originates mainly from red gi- 
ant and supergiant stars, it is present in stars of types A to M. 
Because the li^ht from the galaxy is dominated by the bright 
and young SSCs (I Johnson et al.ll200Qh . it is not the presumably 
relaxed old underlying stellar population that is sampled when 
observing the Ca triplet, but instead fairly young stars — red 
supergiants appear after 5 Myr — that have been formed dur- 
ing the current burst of star formation. This is consistent with 
the fact that we do not see any evidence for absorption in the 
Paschen lines (cf. Fig. [2]). 

The stellar velocity field (panel (c) in Fig. [3]) displays no sig- 
nificant velocity gradient across the field, and thereby no sign 
of rotation. We find the stellar component of He 2-10 to have a 
systemic velocity of 872 ± 6 km s"^ in good agreeme nt with the 
usually adopted systemic velocity for He 2-10 (e.g. JjohanssonI 
11987 ). The absence of a gradient in the stellar velocity field sug- 
gests that the young component, which dominates the light, is re- 
laxed and governed by random motions. The stellar velocity dis- 
persion map (panel (d)) shows an average value of 45 ± 4 km s"^ 
outside the bright starburst centre (region A), while in the very 
centre itself the stellar velocity dispersion drops to 28 ± 3 km s~^ . 
This can be understood by assuming that the young stars, domi- 
nating the light, are being formed in a locally dense and dynami- 
cally less hot environment. Outside the burst region, the stars are 
older and have had more time to virialize with the starburst host, 
thereby increasing the velocity dispersion. 

Assuming virialisation, we can estimate masses using the 
relation Mg- = 1.1 x 10^ • • cr^, where the eff'ective ra- 
dius Tg is given in i inits of kpc, cr in k ms~\ and Mg- in M© 
teender et al.lll992l: iGuzman et aD[l996h . For the central star- 
burst with a size of l.V (derived from archival Hubble Space 



(a) Continuum near CaT (b) [SIM] Emission 




(c) CaT Velocity Field (d) CaT Velocity Dispersion 




(e) [Sill] Velocity Field (f) [Sill] Velocity Dispersion 
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Fig. 3. Panels showing the field of view of ARGUS, 22 by 14 
pixels, corresponding to 1174 x 7! 3. As indicated in Fig.[TJ left 
is north-east and up is north-west. Panels (a) and (b) show the 
continuum emission near the Ca triplet and the monochromatic 
emission from [S iii] respectively, in arbitrary flux units. Panels 
(c) and (d): Heliocentric line-of-sight velocity and velocity dis- 
persion from the Ca triplet. Colours for these and the following 
panels indicate values according to the colour bars beloweach 
column. Panels (e) and (f): Velocity and velocity dispersion from 
[S III]. Panels (g) and (h): Velocity and velocity dispersion from 
the Paschen lines. 



Telescope NICMOS camera //-band photometry) and cr = 28 
kms"^ we calculate the mass to be 4 x 10^ M©. The same exer- 
cise for the region outside region A, represented by an efl'ective 
radius of ^ 6'' and the average stellar velocity dispersion cr = 45 
kms~^ gives a total galaxy mass of 6 x lO^M©. This compares 
very well with the total mass in atomic and molecular gas mea- 
sured by Kobulnicky et al. (1995). The irregular structure of the 
central morphology and possible variations in the mass to light 
ratio leads to an uncertainty in of about a factor of 2, translat- 
ing li nearly to the d erived mass. By integrating the light profiles 
from lNoeske et al.l (12003), we estimate a mass-to-light ratio in 
the H band of 0.05 and 0.2 for the two masses derived above, 
respectively. This is consistent with the very young stellar popu- 
lation in He 2-10. 
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3.2. Gas 

The emission-line velocity dispersions of the [S iii] and Paschen 
emission-lines are almost constant over the field of view, with 
average values 49 + 6 (standard deviation) and 44 + 3 kms"^ 
respectively. In the northern corner of the field, there is an in- 
dication of an increase in cr for the [S iii] and a decrease for 
the Paschen lines, although this is the region with the weak- 
est signal. Overall, these values are higher than the stellar ve- 
locity dispersion at the very centre, but they agree with the 
value outside the starburst region. The observed gaseous veloc- 
ity dispersion is more than a factor two higher than the average 
measured turbul ent velocity dispersion of H ii regions in spiral 
galaxies (e.g., Chu & Kennicuttl ll994') but not unusually high 
for a BC G where it is often used empirically as a mass mea- 
sure fe.g. jMelnick et al.|[T987l: iGuzman eTaDll996l:lQstlin et al.1 
[2QQlh . Since the gas velocity dispersion in He 2-10 is higher than 
that of the stars in the centre, the outflow (see below) probably 
also contributes to the line width. 

The velocity fields of the [S iii] and Paschen emission-lines 
(panels (e) and (g) in Fig. [3]) clearly deviate from that of the 
stars. We observe gas velocities to be generally higher than the 
stellar velocity and find the region west and north-west of the 
central starburst to display velocities reaching up to 50 kms"^ 
higher than the systemic velocity. At the east and south-east, 
the [S III] line-of-sight velocity is less than the systemic value 
by 10 kms"^ Except in regions with barely detected Paschen 
emission, the velocity field of the [S iii] line resembles that of 
the Paschen lines, in agreement with Ha velocities measured by 
iMendez et all (Il999l) in the sense that we observe blueshifted 
lines towards the north-east of the field. W e also find agreement 
with our Ha Fabry-Perot measurements (iMarquart et al.l [20061 
and 2007 in preparation). 

Althou gh on larger scales, the gas in He 2-10 is k nown 
to rotate ( Kobulnicky et aL] 1 19951: iMarquart et"aD [2006), the 
emission-line velocity fields we observe here alone cannot es- 
tablish whether the high and low values to either side of the 
zero-velocity curve are signatures of rotation in the centre of 
He 2-10. The rough orientation of the velocity gradient with 
the morphological major axis might suggest this, but the corre- 
sponding mass from assuming Keplerian motions is an order of 
magnitude below the value derived above. In addition, the lack 
of rotation in the stellar component seems t o preclude the exis- 
tence of a dominant star-forming disk. Kobu lnicky et al.l (1 19951) 
find an elongated cloud of molecular gas with a tail towards the 
south-east (see also our Fig. [T]). Over the region of our field of 
view, the CO gas displays a velocity gradient that is in agreement 
with the gaseous motions presented here. Since the stars do not 
follow these motions, they probably have not formed out of the 
same gas, which might instead represent the part of the cloud 
with high angular momentum or the infall of fresh gas. 

On the other hand, the ISM is subject to strong feedback 
from the starburst that d rives the well- studied outflow in He 2- 
10. Both iBeck & Kov d fl999) and Mendez et al. (1999) have 
observed bubble structures in deep Ha images (cf. Fig.[T]) blow- 
ing o ut towards the n orth-east and south-west of the central star- 
burst. I Johnson et al.l (2000) find this scenario to be consistent 
with the measured energy output from the SSCs in the central 
starburst. We measure blueshifted lines towards the proposed 
outflow regions on both sides that can be attributed to a bipo- 
lar outflow, seen close to edge-on. 

Sch wartz et alj (|2006l) have used ultraviolet (UV) absorption 
lines to measure the speed of the outflow at the position of the 
brightest SSCs and found it to be -170 kms"^ with respect to 



the systemic velocity. In contrast to this, we measure gaseous 
velocities at the systemic value. We attribute this discrepancy 
to the difl'erent regions probed. UV absorption lines sample the 
cold and warm ISM in the foreground, while emission lines trace 
the ionised ISM at the location of the starbrust. 

4. Conclusions 

We have used integral field spectroscopy to measure the stellar 
kinematics in He 2-10 using the Ca triplet absorption lines, and 
compared it with the motions of the ionised gas. We find a large 
discrepancy between stars and gas. The stellar component shows 
no signatures of rotation or any other non-random motions. The 
gas, however, shows bulk motions that are consistent with previ- 
ous results from molecular gas kinematics and the bipolar out- 
flow, driven by the starburst in the centre. 

In particular, IFU spectroscopy has allowed us to show that 
the stars in the centre show no sign of rotation along any direc- 
tion and that gas and stars are kinematically clearly decoupled. 
The virialization of the stellar component is compatible with 
He 2-10 being an advanced merger because simulations show 
that virialization is quick and that the peak in star formation rate 
in a galaxy merging event coincides with the completion of the 
merger (e.g., Barnes & Hernquist 1992; Cox et al. 2006). 

That the stellar component is supported by velocity disper- 
sion rather than rotation also hints toward the future evolution of 
He 2-10 in the direction of a n elliptical- l ike ga laxy (iGeha et al.l 
2002), as earlier suggested by lKoo et"an d 19951) for more distant 
emission line galaxies. Indeed, the strongly peaked light profile 
will probably, after the burst has faded, resemble that of a nucle- 
ated dE like M 32 or NGC 205. The lower velocity dispersion 
that we find for the youngest stellar component in the very cen- 
tre of He 2-10 may in time be dynamically heated and begin to 
resemble the rest of the galaxy. 

Our results also suggest a need for caution in studies of star- 
burst galaxies at higher redshifts, where emission-lines often are 
the only feasible kinematical tracers. Observed spatial variations 
in line-of-sight velocities need not be representative of the grav- 
itational potential if decoupling between gas and stars is com- 
mon. 
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